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Cd , Abstract. We study the inverse spectral problem of reconstruct- 

ing energy-dependent Sturm-Liouville equations from their Dirich- 
let spectra and sequences of the norming constants. For the class of 
problems under consideration, we give a complete description of the 
corresponding spectral data, suggest a reconstruction algorithm, 
and establish uniqueness of reconstruction. The approach is based 
Ph I on connection between spectral problems for energy-dependent 

^^ • Sturm-Liouville equations and for Dirac operators of special form. 

-)— > 

1. Introduction 

The main aim of the paper is to study the inverse spectral problem 
of reconstructing Sturm-Liouville differential equations on (0, 1) with 
l/^ ' energy-dependent potentials from their Dirichlet spectra and suitably 

defined norming constants. The spectral problem of interest is given 
by the differential equation 



OO 

rn 

O' (1.1) -y" + qy + 2\py = \^y 

CN ■ 

and the Dirichlet boundary conditions 

(1.2) y{0) = y{l) = 0. 

Here p is a real-valued function in ^2(0, 1) and g is a real-valued dis- 
c^_' tribution in iy2~^(0! 1)? see detailed definitions in the next section. 

Sturm-Liouville spectral problems with potentials depending on the 
spectral parameter arise in various models of quantum and classical 
mechanics. For instance, to this form can be reduced the correspond- 
ing evolution equations (such as the Klein-Gordon equation [22||39]) 
that are used to model interactions between colliding relativistic spin- 
less particles. Then A^ is related to the energy of the system, thus 
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explaining the term "energy-dependent" in the title of the paper. An- 
other typical example is related to vibrations of mechanical systems in 
viscous media, see [5^ . 

Problems of the form (11. ip have also appeared in the physical litera- 
ture in the context of scattering of waves and particles. In particular, 
M. Jaulent and C. Jean in [T81 - I2T] studied the inverse scattering prob- 
lems for energy- dependent Schrodinger operators on the line; see also 
the papers [IlE3l[iE21ESlESlllIlEIllS2]. An interesting approach to 
the spectral analysis of the Klein-Gordon equations using the Krein 
spaces (i.e., spaces with indefinite scalar products) was suggested by 
P. Jonas 23 and H. Langer, B. Najman, and C. Tretter [SSHMIES]. 

Non-linear dependence of equation (11. II) on the spectral parameter A 
suggests that (ll.ip - (ll.2p should be regarded as a spectral problem for 
a quadratic operator pencil. Although some spectral properties of such 
a problem can easily be derived from the general spectral theory of 
polynomial operator pencils |3l], there have been rather few papers in- 
vestigating the inverse problem of reconstructing the potentials p and q 
from the suitably defined spectral data. The problem with p G W^2^(0, 1) 
and q G i^2(0, 1) and with Robin boundary conditions was discussed 
by M. Gasymov and G. Guseinov in their short paper (TD] of 1981 con- 
taining no proofs. Such problems for (quasi)-periodic boundary con- 
ditions were considered in, e.g., [l2l[T3l|371|38l[55]; however, typically 
only Borg-type uniqueness results were established therein. Some non- 
classical settings of the inverse spectral problem (e.g., those with mixed 
given data, Hochstadt-type problems, or inverse nodal problems) were 
discussed in [2SH2S]. 

The main aim of the present paper is to investigate in detail the 
inverse spectral problem for equations (ll.ip - (ll.2p . under minimal 
smoothness assumptions on real-valued potentials p and q. In par- 
ticular, the distributional potential q can include e. g. the Dirac delta- 
functions of Coulomb-like singularities that are widely used in quantum 
mechanics to model interactions in molecules and atoms, see the mono- 
graphs by S. Albeverio et al. [2] and by S. Albeverio and P. Kurasov [5] 
and the extensive references lists therein. For this wide class of prob- 
lems, we shall give a complete description of the corresponding spectral 
data, suggest a reconstruction algorithm, and establish uniqueness of 
reconstruction. Our approach consists in reducing the spectral problem 
for (ll.ip -( [L2|) to the one for a related Dirac operator; we then study 
possibility of reconstructing a Dirac operator of special form from its 
spectral data. Since various steps of our considerations are explicit, 
they form basis for a reconstruction algorithm. 



The paper is organized as follows. In the next section, we introduce 
the main objects of study and formulate the main results. In Section [31 
we show that the spectral problem fll.ip - fll.2p is closely related to the 
one for a Dirac operator of a special form. Further we discuss the pos- 
sibility to transform Dirac operators to some canonical form by means 
of the so-called transformation operators. Existence and uniqueness 
of a special Dirac operator (and thus of the corresponding operator 
pencil) possessing the prescribed spectral data is tackled with in Sec- 
tions and E] respectively. Finally, the last Section [7] summarizes the 
reconstruction algorithm and discusses some possible extensions. 

Notations. Throughout the paper, L2,r(0, 1) and Vr2~]^(0, 1) will 
stand for the sets of real- valued functions in 1/2(0, 1) and distributions 
in W2^{Q, 1), respectively. We denote by p{T) and cr(T) the resolvent 
set and the spectrum of a linear operator or a quadratic operator pen- 
cil T, and by A^2 = A^2(C) the linear space of 2 x 2 matrices with com- 
plex entries endowed with the Euclidean operator norm. The super- 
script t designates transposition of vectors and matrices, e.g., (ci,C2)* 
is the column vector (^^) . 

2. Preliminaries and main results 

Equation (11.11) contains terms depending both on A and A^ and there- 
fore leads to a spectral problem for a quadratic operator pencil that we 
now introduce. 

To begin with, we recall that the potential q in (II. ID is a real- 
valued distribution in the Sobolev space W2^{Q,1) and thus q = r' 
for some r G L2,ir(0, 1). The Sturm-Liouville operator with potential q 
can be defined following the regularization method due to Savchuk and 
Shkalikov [l5lll6]. Namely, for every absolutely continuous function y, 
we denote by y'-^' := y' — ry its quasi-derivative and introduce the 
differential expression 

m:=-{y^'^y-ry^'^-r'y 

acting on 

domi ={ye AC{0, 1) | yW E AC[0, 1], iiy) E 1^(0, 1)}. 

We now define the operator A via 

Ay = iiy) 

on the domain 

domA := {y E dom£ | y{0) = y{l) = 0}. 
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A straightforward verification shows that £(y) = —y" + qy in the 
sense of distributions and, if q is integrable, even in the usual sense. 
Therefore for regular g, A is the standard Sturm-Liouville operator 
with potential q and Dirichlet boundary conditions. It is known [^B6] 
that if g G W2^{^, 1) is real- valued, then the operator A is self-adjoint, 
bounded below, and has a simple discrete spectrum. The Green func- 
tion of the operator A is continuous in the square [0, 1] x [0, 1], so that 
the resolvent of A is of Hilbert-Schmidt class. 

Further, we denote by B the operator of multiplication by the po- 
tential p G 1^2(0, 1). The operator B is in general unbounded; however, 
since BA~^ is of Hilbert-Schmidt class, B is A-compact [25, Ch. IV]. 
In particular, dom B D dom A and B is bounded relative to A with rel- 
ative A-bound [HI Lemma III. 2. 16]. Finally, / stands for the identity 
operator in ^2(0, 1). 

Now the spectral problem f ll.ip - fll.2p can be regarded as the spectral 
problem for the quadratic operator pencil Tp^q defined as 

Tp,q{X) := X^I - 2XB - A 

for A G C on the A-independent domain doniTp^ := dom A. The 
above-mentioned properties of the operators A and B imply that, for 
every A G C, the operator Tp^q{\) is well defined and closed on dom Tp^g. 
Before continuing, we recall the following notions of the spectral theory 
of operator pencils, see |3^ 



Definition 2.1. The spectrum a{Tp^q) of the operator pencil Tp^q is the 
set of all A G C for which Tp_g(A) is not invertible, i.e., 

CT{Tp,q) = {A G C I G CT{Tp,q{X))}. 

A number A G C is called the eigenvalue of Tp^q if Tp^q{X)y = for 
some non-zero function y G domTp^, which is then the corresponding 
eigenf unction. Finally, 

p{Tp^q) := C \ a{Tp^q) 

is the resolvent set of the operator pencil Tp^q. 

It was shown in [42] that the spectrum of the pencil Tpq consists 
entirely of eigenvalues and that criTp^q) is a discrete subset of C. In 
general, Tp ^ can possess non-real and/or non-simple eigenvalues; the 
latter means that the algebraic multiplicity of an eigenvalue can be 
greater than 1 [33]. The approach to the inverse spectral problem 
for Tp^q we are going to use is currently well understood only for the 
case where the spectrum of Tp^q is real and simple. This happens e. g. 
when there is a Aq G M such that the operator Tp_q(Ao) is negative [42] . 
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i.e., when Tp^g is a hyperbolic pencil [Ml Ch. 31]. Making the shift 
A I— 7- A + Ao if necessary, we may (and shall) assume that Aq = 0. 
Therefore, our standing assumption throughout the paper is that 

(A) p G L2,k(0, 1), q G W^2"r(0' ■^); ^^^ ^^^ operator A is positive. 

Under assumption (A), the eigenvalues of Tpg are all real, simple, 
and can be labeled in increasing order as A„ for n G Z* := Z \ {0} so 
that 

A„ = 7m + po + \n, 

where po '■= /g p{^) dx and (A„) is a sequence in £2(Z*) [lUP^] . Observe 
also that under assumption (A) the point A = is in the resolvent set 

of Tp^q. 

We observe that equation f 1 1.11) can be recast using quasi-derivatives 
as 

(.{y) + 2\py = X^y 

and that for every complex a and b it possesses a unique solution sat- 
isfying the initial conditions y{0) = a and yf^'(O) = b. This allows to 
introduce the norming constants in the following way. 

Definition 2.2. For an eigenvalue A„ of Tp^g, denote by yn the corre- 
sponding eigenfunction normalized by the initial conditions yn{0) = 
and yii (0) = A„. Then the quantity 

(2.1) an -.= 2 I yl{t)dt - |- / p{t)yl{t)dt 

Jo ^n Jo 

is called the norming constant corresponding to the eigenvalue A„. 

Although this definition looks somewhat artificial, there is a good 
reason for defining the norming constants via (12. ip . Firstly, for p = 
0, problem (ll.ip -( !L2|) becomes the spectral problem for the Sturm- 
Liouville operator A, and (12. ip up to a constant factor agrees with the 
common definition of the norming constant [llj. Secondly, denoting by 
Tpg the A-derivative of Tp^g, we find that 

whence a„ determines the type of the eigenvalue A^, see [31]. In the next 
section, we shall transform the spectral problem for the operator pencil 
to that for some Dirac operator ^{P), and (12. ip agrees with the stan- 
dard definition of the norming constant for ^(P). Next, we mention 
that the function u{x, t) := ?/„(x)e^"* is a solution of the corresponding 
evolution equation Tp^q{d/dt)u = 0, and its energy {Au, u) + {u, ii), with 
it denoting du/dt, is equal to A^a^e^'*'"*. 
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Remark 2.3. It follows from |12] that if A is invertible, then its posi- 
tivity (cf. assumption (A)) is necessary and sufficient for all a^ to be 
positive. 

Definition 2.4. Assume that A is an eigenvalue of the quadratic op- 
erator pencil Tp^q and that a is the corresponding norming constant. 
Then (A, a) is called the spectral eigenpair of Tp^q. The spectral data 
sd(Tp_g) of the pencil Tp g is the set 

sd(TpJ:={(A,a)| AG(7(TpJ} 

of all its spectral eigenpairs. 

The inverse spectral problem of interest is to reconstruct the poten- 
tials p and q of the operator pencil Tp^q given its spectral data sd(Tp g). 
Some properties of spectral data for the class of quadratic pencils Tp^q 
under consideration was established in [12]. Our aim in this paper is, 
firstly, to give a complete description of the set of the spectral data and, 
secondly, to find and justify an algorithm reconstructing the potentials 
p and q from the spectral data. 

We observe that when p = 0, the spectral problem for the operator 
pencil To^q becomes the spectral problem Ay = \^y for the Sturm- 
Liouville operator A. Then A_„ = — A„, «_« = Oi^ and «„ of (12. ip 
agrees with the standard definition of a norming constant [TT]. For the 
Sturm-Liouville operator A with a real- valued q G i^i(0, 1) and Robin 
boundary conditions, it was proved in [TT] that the spectrum (A^)jigN 
of A and the sequence (a„)nGN of the corresponding norming constants 
uniquely determine the potential g; the case of a distributional poten- 
tial q G W2^{Q, 1) was treated in e.g. |6l[T6lll7]. The above references 
also suggest algorithms of reconstructing the potential q from the spec- 
tral data of A. 

The operator pencil Tpq contains two real-valued potentials p and 
q to be determined in the inverse problem; however, the spectral data 
of Tp^q are twice as large as for a standard Sturm-Liouville operator. 
Therefore one may hope that the inverse spectral problem of recon- 
structing p and q from the spectral data of Tpq is well posed. Our first 
main result gives uniqueness of reconstruction. 

Theorem 2.5. Under assumption (A), the operator pencil Tp^q is 
uniquely determined by its spectral data sd(Tp_g). 

It follows from ^^ (cf. also the result of [lO] for p G IVsHO, 1) 
and q G L2(0, 1)) that the spectral data for the operator pencils under 
consideration belong to the following set. 



Definition 2.6. We denote by SD the family of all sets {(A„, a„)}„gz* 
consisting of pairs (A„, a„) of real numbers satisfying the following 
properties: 

(i) A„ are nonzero, strictly increase with n G Z*, and have the 
representation A„ = vm + /i + A„ for some /i G M and a sequence 
(A„) in £2(Z*); 
(ii) a„ > for all n ^ Z* and the numbers a„ := «„ — 1 form an 
^2(Z*)-sequence. 

Our second result claims that, conversely, every element of SD is 
spectral data for some operator pencils Tp^g under consideration. In par- 
ticular, it shows that conditions (i)-(ii) above give a complete descrip- 
tion of the spectral data for the operator pencils Tp^g with p G L2,r(0, 1) 
and gG 1^2:^(0,1). 

Theorem 2.7. For every sd G SD, there exists p G L2,k(0, 1) and q G 

V1^2~r(05 1) such that sd is the spectral data for the operator pencil Tp^g. 

The proof of this theorem is constructive and suggests the explicit 
reconstruction algorithm determining the potentials p and q from the 
set sd belonging to SD; see Section [71 

Our approach consists in reducing the spectral problem for Tp^g to the 
one for a Dirac operator of a special form acting in L2(0, 1) x L2(0, 1), 
see Section [21 Under a suitable unitary gauge transformation this Dirac 
operator takes the "shifted" AKNS normal form. For AKNS Dirac 
operators, the direct and inverse spectral problems are well understood, 
see [3l[9l[3T]. We shall use the known methods to first reconstruct 
the Dirac operator in the "shifted" AKNS form from the given data 
and then to transform this Dirac operator to the one that is directly 
associated with some pencil Tp^g. The latter gives the required solution 
of the inverse spectral problem of interest. 

3. Reduction to the Dirac system 

In this section we shall show that under the standing assumption (A) 
the spectral problem for the operator pencil Tp ^ can be reduced to the 
one for a special Dirac operator. We start with the following observa- 
tion. 

Lemma 3.1. Under the standing assumption (A) the equation i{y) = 
possesses a solution y that is strictly positive on [0, 1]. 

Proof. For /i G M, consider an operator A^ defined via A^y = £{y) on 
the set of functions in dom i verifying the boundary conditions 

y(0)-/iyW(0)=y(l)=0. 



The operators Ah are self-adjoint and have discrete spectra. We show 
that Ah depend continuously on h in the norm resolvent sense. 

Indeed, fix a non-real A and denote by (/?_ and ip^ solutions of the 
equation i{y) = \y satisfying the initial conditions v'_(0) = 0, V'L (0) = 

1 and the terminal conditions (p+{l) = 0, V'tj. (1) = 1, respectively. Take 
an arbitrary / G -^2(0, 1) and set gh '■= {Ah — A)^^/. The function 
z := Qh — go solves the equation i{z) = Xz and satisfies the terminal 
condition z{l) = 0; therefore there exists a constant Ch{f) such that 
-2 = Ch{f)<p+- Since the function Qh = Qo + Ch{f)'^+ satisfies the initial 
condition ghiO) = hg\^ (0), we conclude that 

e,(/) = "go"'") 

"'" ^.(0)-fc^!;l(0)- 

The denominator never vanishes for /i G M since the non-real A is in 
the resolvent set of all operators Ah- 

With W := v^+v'L — V'+V- denoting the modified Wronskian, we 
find that W is constant on [0, 1] and that 

9o{x) = ^ £ V4t)m dt + ^ £ V4t)f{t) dt, 
whence 

9ll\o) = ^J\4t)f{t)dt. 

Therefore the resolvent difference {Ah — A)~^ — {Aq — A)~^ is a rank 
one operator equal to 



and thus the norm resolvent continuity of Ah as a function of /i G M 
follows. 

By assumption (A), the operator A = Aq is uniformly positive; thus 
the above continuity in h implies that there is an /i > such that Ah ^ 
0. Fix one such an h and denote by i/h the solution of the equation 
i{y) = satisfying the initial conditions y{0) = h and y^^^{0) = 1. We 
claim that yh stays positive over [0, 1]. 

To prove this, we consider the quadratic form ah of the operator Ah. 
Integration by parts shows that 

ah[y] = \\y'\\^-2Reiry',y) + ^\y{0)\-' 
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for y G dom Afi- Next (cf. [15]), the quadratic form Re{ry',y) is rela- 
tively bounded with respect to the form 

My] ■■= bt + l\ym' 

with relative bound 0. Since the quadratic form O-h is closed on the 
domain 

dom~aH:={yeW^{0,l)\y{l) = 0}, 

Theorem VI. 1.33 of [25j implies that the quadratic form ah is closed on 
the same domain. 

Now assume, on the contrary, that the above function yh does not 
remain positive over [0, 1] and denote by xq G (0, 1] a zero of yt- Then 
the function 

\yh{x) for X < xq] 
z{x) := < 

10 for X > Xo 

belongs to the domain of the quadratic form ah of Ah- Moreover, 
integration by parts in the expression 



\y'h{x)\^ dx - 2Re / {ry'jJJh){x) dx + -\yh{Q)\ 
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shows that ah[z\ = 0, which contradicts uniform positivity of the op- 
erator Ah- The derived contradiction shows the assumption on yh was 
wrong, and thus yh is a solution of i{y) = that is positive over [0, 1]. 
The proof is complete. D 

Now we take a solution y of the equation i{y) = that stays positive 
on [0, 1] and set v := y'/y. Since y' = y^^^+ry G i^2(0, 1), the function v 
is in -^2(0, 1). Moreover, direct calculations show that q = v' + v"^ (so 
that g is a Miura potential, see [21]) and that the operator A can be 
written in the factorized form, viz. 

P'i> ^^ = -(5^ + ")(^ -")■''■ 

We observe that there are many different v satisfying the Riccati equa- 
tion v' + v'^ = q and thus allowing the above factorization of A; we 
fix one such a f in what follows and notice that the function v — r is 
continuous. 

For A 7^ we can recast the spectral problem (11.11) as a first order 
system for the functions U2 '■= y and Ui := {y' — vy)/X, namely, 

(3.2) U2 — VU2 = Xui, 

(3.3) —u[ — vui + 2pu2 = Xu2. 

9 



Setting 

we see that the above system is the spectral problem ^{P)u = Au for 
the Dirac operator ^(P) acting in L2(0, 1) x L2(0, 1) via 

(3.5) ^(P)u = J^ + Pu =: £(P)u 

ax 

on the domain 

dom^(P) :={u= {ui,U2y eW^{0,l)xW^{0,l) | ^2(0) = M2(1) = 0}. 

Lemma 3.2. The nonzero spectra of the Dirac operator ^(P) and the 
operator pencil Tp^g coincide. 

Proof. We start by observing that the Dirac operator &{P) is self- 
adjoint in the Hilbert space Iv2(0, 1) x 1/2(0, 1) and has a simple discrete 
spectrum [21]. The above arguments show that every eigenvalue A 
of Tp^q is an eigenvalue of ^(P) as well. 

Conversely, assume that A 7^ is an eigenvalue of S>{P) and 
u = (mi, U2)* is a corresponding eigenfunction. Determining ui via 
U2 from (13. 2 p and substituting in (13. 3p . we find that y := ti2 satisfies 
the equality 

which in view of (13.11) yields Tp^qH = 0. Clearly, y is non-trivial and 
satisfies the Dirichlet boundary conditions; henceforth it is an eigen- 
function of Tp q corresponding to the eigenvalue A. D 

Remark 3.3. A straightforward analysis of (I3.2p -( 1X^ shows that A = 
is an eigenvalue of the Dirac operator ^(P), the corresponding eigen- 
function being u = (ui,M2)* with ui = exp{— J v) and M2 = 0. How- 
ever, under assumption (A) the number A = is never an eigenvalue 
of Tp q-, therefore, 

a(^(P)) = a{Tp,q) U {0}. 

The norming constant for the Dirac operator ^(P) corresponding to 
an eigenvalue A is defined as 

iiuir = ikiiii^ + \\u2\\i„ 

where u = {ui,U2Y is the eigenfunction for A normalized by the initial 
conditions mi(0) = 1 and ^2(0) = 0. Assume that A 7^ 0; then, as 
shown in the proof of Lemma 13. 2[ y := U2 is an eigenfunction of Tp g 
corresponding to the eigenvalue A; note also that y is real valued and 
satisfies the initial conditions y{0) = and yt^'(O) = (y' — ry)(0) = 
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iy' — vy){0) = Ami(0) = A. Integration by parts on account of (13.11) 
gives 

lk2-^W2||i2 = {Au2,U2)l2; 

using now ( 13.2^ . we conclude that 

||uf = —{Ay, y) + {y, y) = -j{By, y) + 2{y, y), 

which coincides with f l2.ip . We have thus estabhshed the following 
important result. 

Lemma 3.4. The norming constants corresponding to nonzero eigen- 
values of the Dirac operator &{P) and the operator pencil Tp^g coincide. 

The above lemma suggests that we can use the spectral data sd(Tp^q) 
of the operator pencil Tp^q in order to find the related Dirac opera- 
tor ^(P). Having determined the potential P = {pij)'ij=i of ^(-P), 
we then identify the potentials p and q of the operator pencil Tp^g as 
P '■= P22/2 and q := —p'12 + Pu- We note, however, that since the 
factorization (13.1 p is not unique, there are many Dirac operators &{P) 
associated with Tp^g. Therefore the spectral data sd(Tp,j) cannot de- 
termine such an operator ^{P) uniquely. 

The reason for this non-uniqueness is quite clear from Remark 



and Lemma 13. 4t indeed, the spectral data for Tp^g leave the norming 
constant ao for the eigenvalue Aq := of ^(P) undetermined. It is this 
freedom in the choice of ao that leads to non-uniqueness of potentials P 
for the associated Dirac operators ^(P). However, we shall show in 
Section [H] that all such Dirac operators determine the same pencil Tp^g. 

4. Transformation operators 

We shall use the relation between the operator pencil Tp^g and the 
Dirac operator ^(P) of fl3.4p - fl3.5p explained in the previous section 
and the well-developed inverse spectral theory for Dirac operators to 
reconstruct ^(P) from the given spectral data. Once such a Dirac 
operator has been found, it is then straightforward to determine the 
corresponding potentials p and q of the pencil Tp g. 

However, the classical inverse spectral theory reconstructs a Dirac 
operator with potential in the AKNS form or in other canonical form, 
see Section [51 Therefore we then have to transform such a canonical 
Dirac operator to the one of the form (13.40 keeping the spectral data 
unchanged. This is done by means of the so-called transformation 
operators, which we study in this section. 
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More exactly, assume that P and Q are 2x2 matrix- valued potentials 
in L2((0, 1), A^2) and set 

V^ ■■= {{ui, u^y e W^iO, 1) X W^iO, 1) I ^2(0) = 0}. 

We need a transformation operator ^ = J^{P,Q) between the Dirac 
operators i{P) and £{Q) acting on the set Vq, i.e., for an operator 
satisfying the relation ^i{P)u = i{Q)^u for all u G Pq- Similar 
transformation operators for P and Q with Lipschitz continuous en- 
tries were constructed in [3131]; it is thus reasonable to look for the 
transformation operator <^ of a similar form 



(4.1) ^u{x) = R{x)u{x)+ K{x,s)u{s)ds, 

Jo 

where R and K are 2x2 matrix-valued functions of one and two 

variables respectively. Keeping in mind that the Dirac operators of 

interest, i{P) and i{Q), are considered on functions satisfying the same 

initial conditions, we impose the restriction -R(O) = / guaranteeing 

that ^ preserves the values of functions at x = 0. Under such a 

normalization, R will explicitly be given as 

(42) R{x) = e^'^''U ^°^^2(a;) s\ne2{x) \ ^ ^e,{x)i+e2(x).j 
^ ' ' ^ ' V — sin6'2(x) cos6'2(x) / ' 



with 



(4.3) 



Oi{x) = ^j tT[J{Q{s)-P{s))]ds, 
92{x) = ^j\T{Q{s)~P{s))ds, 



cf. [7]. More exactly, the following analogue of Theorem 3.1 of [7] holds 

true. 

Theorem 4.1. Assume that P and Q are in L2((0, 1), A^2)- Then 
an operator J^{P,Q) of the form (14. ip . with R obeying the condi- 
tion -R(O) = / and with a summable kernel K , is a transformation 
operator for i{P) and i{Q) on the set Vq if and only if the matrix- 
valued function R is given by (14.21) - (14. 3p and the kernel K = {Kij)fj^-^ 
is a mild solution of the partial differential equation 

(4.4) Jd,K{x, y) + dyK{x, y)J = K{x, y)P{y) - Q{x)K{x, y) 

in the domain Vt = {{x,y) \0<y<x<l} satisfying for < x < 1 
the boundary conditions 

(4.5) K{x, x)J - JK{x, x) = JR'{x) + Q{x)R(x) - R{x)P{x), 

(4.6) Ku{x,0) = K22{x,0) = 0. 
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The proof of this theorem follows in general that of Theorem 3.1 
of [7]. One essential difference is that since the matrix- valued func- 
tions P and Q are less regular, K need not be differentiable in the 
usual sense. Therefore differentiation of K should be understood in 
the distributional sense; that is why K is only required to be a mild 
solution of (14.41) . For the sake of completeness, we justify the steps 
involving differentiation in the proof below. 

Proof. Firstly, we observe that, for an integrable K and for u G Pqi 
the standard formula 

— / K{x,s)u{s)ds = K{x,x)u{x) + / dxK{x,s)u{s)ds 
dx Jq Jq 

remain to hold in the sense of distributions, so that 
(4.7) 
^{Q)3^u{x) = JR{x)u'{x) + {JR\x) + Q{x)R{x)}u{x) 

+ J-f- / K{x,s)u{s)ds + / Q{x)K{x,s)u{s)ds 
dx Jo Jo 

= JR{x)u'{x) + {JR'{x) + Q{x)R{x) + JK{x,x)}u{x) 



+ / {Jd^K{x,s) + Q{x)K{x,s)}u{s)ds. 
Jo 

Similarly, the integration by parts formula 

PX PX 

/ K{x,s)u'{s) ds = K{x,x)u{x)—K{x,0)u{0)— dsK{x, s)u{s) ds 
Jo Jo 

holds in the distributional sense, and in the same sense one gets the 
equalities 

(4.8) 

^i{P)u{x) = R{x)Ju'{x) + R{x)P{x)u{x) 

+ f K{x,s){Ju'{s) + P{s)u{s)}ds 
Jo 

= R{x)Ju'{x) + {R{x)P{x) + K{x, x)J}u{x) 

PX 

-K{x,0)Ju{0)+ / {K{x,s)P{s)~dsK{x,s)J}u{s)ds. 



JO 

If R is given by (giS]) and K satisfies dH]), then JR = RJ and 
K{x,0)Ju{0) = for u G Vq; using now (14.41) and (14. 5 P in equations 
(14.71) and (14. 8p . we arrive at the equality 

e{Q)^u= jri{P)u 

on the domain Vq, thus establishing the sufficiency part. 
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To prove necessity, we equate the integrands and the coefficients of u 
and u' in fl4.7p and ( 14. 8p . The coefficients of u' yield the relation JR = 
RJ and so 

a(x) b{x) 



^ —b{x) a{x) 

Now equate the coefficients of u to get the relation 

(4.9) JR\x) + Q{x)R{x) - R{x)P{x) = K{x, x)J - JK{x, x) 

coinciding with (14. 5p . Using the fact that KJ — JK and JKJ + K 
have zero traces, we derive from (14. 9 p the following system for a and h: 

aim d(a\_(ivJ{Q-P) -tTiQ-P)\fa 

^ ^^ dx[b )-[ tr(Q-P) tr J(Q-P) ) [ b 

Multiplying the ffist row by a and the second by b and adding gives 

{a^ + b^y = {a' + b^)tTj{Q-P), 

so that a^ + 6^ = cexp(26i) with 6i of ( 14.31) and some constant c. 
Recalling the assumption -R(O) = /, we conclude that c = 1. 

Now, upon substituting a = exp(6'i) cos?7 and b = exp(6'i) sinr^ in the 
system (I4.10p . we find that ?7 = 6*2 + Ci with a constant Ci. Using again 
the normalization -R(O) = J, we obtain rj = 62 + 27m, n G Z. Thus the 
matrix R is indeed given by (14. 2p . 

Next, the equation for the kernel K follows from the equality of the 
integrands in (14. 7p and (14. 8p . Finally, as the term i^(x, 0) Ju(0) must 
vanish for all u G Vq, relation (14.61) follows. The proof is complete. D 

The above theorem reduces the question on existence of the trans- 
formation operator to that on solvability of the system (I4.4l) - (l4.6p . 
Existence of mild solutions to that system is discussed in the next the- 
orem. 

Theorem 4.2. Assume that matrix-valued functions P and Q are 
in iv2((0, 1), A^2)- Then the system (14. 4p - (14.61) has a unique solution in 
the sense of distributions; moreover, this solution belongs to L2(fi, A^2)- 

Proof. The proof of this theorem is rather standard and uses reduction 
to the equivalent system of integral equations. We only sketch the main 
idea and refer the reader to the paper [53] where the missing details 
can be found. 

Introduce a four- component vector-function L = (Li, L2, L^, L4)* via 

Li = K21 — K12, L2 = K22 — Kn, 
L3 = Kn + K22, L4 = K12 + K21. 
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In terms of this vector the system (14. 4p takes the form 

(4.11) (9,. + Edy)L{x, y) = F{x, y)L{x, y), 

where E = diag(l, — 1, 1, — 1) and F{x,y) is a 4 x 4 matrix-function 
whose entries are hnear combinations of the entries of P{y) and Q{x). 
The boundary conditions (14. 5 p and (14. 6 p for K translate into the rela- 
tions 

Lk{x,x) = gk{x), fc = 2,4, 

(4.12) Li{x,0) = Li{x,0) 

L3(x,0) = -L2(x,0), 

where g2 and g4 are respectively the (1, 2)- and (1, l)-entries of the 
matrix R{x)P{x) — Q{x)R{x) — JR'{x). Under the assumptions of the 
theorem, g2 and g^ belong to -^2(0, 1). 

We next denote by Fj, i = 1,4, the i-th row of the matrix F and 
rewrite the system (l4.1ip -( P:.12p as a system of the integral equations 

Li= Fi{-s + x + y, s)L{-s + x + y, s)ds + gi (— ^j > ^ = 2,4 

ry 
Li = I Fi{s + X — y, s)L{s + X — y, s)ds 
Jo 

x — y , - 

/2 /it ?/ \ 

F4(-s + x-y, s)L{-s + x-y, s)ds + g^ ( j 

F^{s + X — y, s)L{s + X — y, s)ds 

2 / X — V 

F2{-s + x -y, s)L{-s + x -y, s)ds - g2 I 
'0 V ^ 

Applying to the latter the successive approximation method and using 
the fact that (?2 and 5(4 belong to ^2(0,1), one proves existence of a 
unique solution L belonging to L2 (fi, C^) . This gives a unique mild 
solution K of the original hyperbolic system (I4.4p - (l4.6p and proves 
that K eL2{n,M2). □ 

Throughout the rest of the paper, we shall assume that the matrix- 
valued potentials P and Q are Hermitian, i.e., that P*{x) = P{x) 
and Q*{x) = Q{x) a.e. on [0, 1]. Then the corresponding Dirac opera- 
tors ^(P) and ^{Q) are self-adjoint and have simple discrete spectra. 
Moreover, the eigenvalues of f^n{.P) of ^(P) can be labeled by n G Z 
so that Hn = T^n + ^ti P + o(l) as \n\ — >■ 00 (see [31]), and similarly for 
the eigenvalues of ^{Q). 
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Just as for an operator pencil Tp g, we define the spectral data for a 
Dirac operator ^ as tlie set {(A, a) | AG cr{^)} of all eigenpairs (A, a) 
composed of eigenvalues A and the corresponding norming constants a. 
We denote by sd(P) (resp. by sd(Q)) the spectral data of ^(P) (resp., 
the spectral data of ^{Q)). 

Also, ^ = ^{P,Q) will stand for the transformation operator of 
the form (14.11) for the differential expressions i{P) and i{Q) on the 
domain Vq, with R given by (14.21) and fl4.3p . We shall write .^ = 
^ + J^, where ^u{x) := R{x)u{x) is the operator of multiplication 
by R and 



J(fu{x) := / K(x, s)u(s) ds 
Jo 
is the corresponding integral operator. We observe that for Hermit- 
ian P and Q the functions i6i and 62 are real valued; in particular, the 
operator M is unitary. 

Before discussing further properties of the transformation opera- 
tor e^, the following simple but useful remark seems in place. 

Remark 4.3. Since the transformation operator J^ intertwines the 
Dirac differential expressions i{P) and i{Q), it is straightforward to 
see that the relation 

(£(P)-A)u = f 

holds if and only if for v := J^Tu and g := J^Tf one gets 

(^(Q)-A)v = g. 

Lemma 4.4. Assume that P and Q are Hermitian and that sd(P) = 
sd{Q). Then the integral operator ^ in the transformation opera- 
tor ^ = 3^{P,Q) is the zero operator. 

Proof. We enumerate the eigenpairs in sd(P) = sd{Q) as (A„, a„) for 
n & 'L and denote by u„ the eigenfunction of the operator ^(P) cor- 
responding to the eigenvalue A„ and normalized via u„(0) = (1,0)*. 
Then v„ := ^(P, Q)u„ is the corresponding eigenfunction for the op- 
erator ^(Q) satisfying the same initial condition. 

The sequences {yin)n&i and (v„)„gz form orthogonal bases of the 
Hilbert space L2((0, 1), C^); moreover, ||u„|| = ||v„|| = y/a^. Thus we 
conclude that the operator ^ is unitary, i.e., that 



+ ^)*(^ + jr) = ^, 

where J^ is the identity operator in -L2((0, 1), C^). Since M is unitary, 
the last equality may be rewritten as 

(^ + ^~^^)*(^ + ^-\J^) = J. 
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We recall that J^ (and thus =^^^J^) is an integral operators with 
lower-triangular kernel belonging to L2(fi,A^2)- Slightly modifying 
the arguments of [ISj Sec. IV. 1], one sees that ^~^.J^ is a Volterra 
operator, whence the inverse (j^ + ^~^J^)~^ exists and is given by 
the Neumann series, 

(^ + ^-^JT)-^ = J- ^-^X + {^-^Jiff + . . . = ^ + ^ 

with J^ being an integral operator with lower-triangular kernel. On 
the other hand, the operator (^~^J^)* is an integral operator with 
upper-triangular kernel, and the relations 

imply that {JM~^J^Y = J^ = 0. Therefore J^ = 0, and the proof is 
complete. D 

The following theorem gives necessary and sufficient conditions on 
transformation operators ^{P,Q) in order that the Dirac operators 
Si{P) and S'{Q) should have the same spectral data. It will essentially 
be used in the reconstruction procedure of the next section. 

Theorem 4.5. Assume that the matrix potentials P and Q are Her- 
mitian. Then the spectral data for the operators &{P) and ^{Q) co- 
incide, i.e. sd(P) = sd{Q), if and only if the transformation opera- 
tor ^{P,Q) for i{P) and i{Q) on the domain Vq only contains the 
unitary part M (i.e., J^ = 0) and ^2(1) = ^'^ for some n G Z. 

Proof. Necessity. If the spectral data for the operators ^(P) and ^{Q) 
coincide, then J^ = by Lemma [4.41 and thus ^{P,Q) = =^. Take 
an arbitrary eigenvalue A of ^(P) and denote by u = (ui,M2)* the 
corresponding eigenfunction. Then 

/ N ^ (a( \\( cos6'2(x)mi(x) + sin6'2(x)M2(x) 

v(x) := ^u = exp(6'i(x)) . ^\ \ ^,\ , ^/. W 

^ ' ^v \ // \^ — sm6'2(a;)-ui(x) + COS ^2 (3^)^2(3;) 

is the eigenfunction for the operator &{Q) corresponding to the same 
eigenvalue A. Observe now that the second components of u(l) 
and v(l) are equal to zero. Since 

(.-. la l^\\ ( cos6'2(1)mi(1) \ 

^W=^^P(^^(^))^-sin^2(lK(l) J 

and U\{V) 7^ 0, we conclude that sin6'2(l) = 0, and thus that 6'2(1) = vm 
for an integer n as claimed. This completes the proof of the necessity 
parts. 

Sufficiency. Suppose that <^(P, Q) = M and that A is an eigenvalue 
of the operator S>{P) with a corresponding eigenfunction u. Consider 
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the vector v = Mu; then i{Q)v = Av. The assumption 6*2(1) = Tin for 
an integer n imphes that -R(l) is a multiple of the identity matrix /. 
Since also -R(O) = /, we conclude that the second component V2 of v 
vanishes at both endpoints and thus v G doTa!^{Q). 

The above implies that the spectrum of S>{P) is contained in the 
spectrum of &{Q). As M is boundedly invertible, the roles of P and 
Q can be interchanged, thus showing that the spectra of the operators 
coincide. Finally, since the operator M is unitary and preserves the 
initial conditions, the norming constants for Si{P) and &{Q) are equal, 
i.e., sd(P) = sd{Q). The proof is complete. D 

5. Reconstruction of the pencil: Existence 

Our aim in this section is to prove Theorem 12.71 on existence of a 
pencil Tp^q with potentials p E L2,r(0, 1) and q G Vr2"]^(0, 1) having the 
prescribed element sd of SD as its spectral data. The uniqueness of 
reconstruction will be dealt with in the next section. 



Outline of the proof of Theorem 2.7 First, we construct a Dirac op- 



erator &{Q) in the "shifted" AKNS form whose nonzero spectrum 
and corresponding norming constants are given by sd. Then we use 
the transformation operator technique to transform S){Q) to another 
Dirac operator !^{P) with the same spectral data and with poten- 
tial P = (pij) of the form f l3.4p . Setting then 

(5-1) ^■"^' ^-=-^12+^12 

and recalling the results of Section |3l we conclude that the operator 
pencil Tp^q is a solution of the inverse spectral problem, thus completing 
the proof. D 

The rest of this section contains details of constructing the poten- 
tial Q of the "shifted" AKNS form and transforming it to a P of the 
form dSa]). 

We start with taking an arbitrary set sd = {(A„, Q;„)}„gz. in SD. 
The enumeration of A^ is uniquely determined by the requirement that 
A_i < and Ai > and fixes the number h in the asymptotic represen- 
tation of part (i) of Definition 12.61 We then take an arbitrary ao > 0, 
put Ao := 0, and augment the set sd with (Aq, ao) to become sd*. 

Now we recall the following facts from the inverse spectral theory 
for AKNS Dirac operators, see [21EI]- Denote by Qq the set of 2 x 2 
matrix- valued functions Qq of the AKNS normal form, namely 

(5.2) Qo := |Qo = (ll ^l) I Qj G ^2,r(0, 1)| . 
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For Qq G Qo) the Dirac operator S>{Qq) is self-adjoint, has a simple 
discrete spectrum, and its eigenvalues can be enumerated as A„(Qo)) 
n G Z, so that XniQo) increase with n and XniQo) = Ten + \n{Qo), 
with an £2(Z)-sequence (A„). The corresponding norming constants 
«n(Qo) are positive and the remainders <5n(Qo) := «n(<5o) — 1 form an 
£2(Z)-sequence. 

Conversely, it follows from the results of [21I1S] that every set 
{{Xn,C(n)}n€Z with A„ and q;„ possessing the above properties is the 
set of spectral data for a unique AKNS Dirac operator ^(Qo) with 

Qo e Qo- 

The set sd* (i.e., the set sd augmented with the pair (Ao,ao) as 
above) looks like spectral data for an AKNS Dirac operator, save that 
the asymptotics of A„ is shifted by some number h; cf. the definition 
of the set SD. For /i G M, we denote by 

Qh := {Qo + hI\Qoe Qo} 

the set of h- shifted AKNS matrix potentials; then the following holds 
true. 

Proposition 5.1. For an arbitrary sd G SD, fix h E M in the rep- 
resentation of (i) in Definition \2.6\ so that A_i < and Ai > 0, and 
denote by sd* augmentation of sd by a pair (Ao,«o); with Aq := and 
a positive ao. Then there exists a unique potential Q & Qh such that 
sd* gives the spectral data for the Dirac operator !^{Q). 

For an arbitrary Q G 1^2 ((0, 1), A^2)5 we introduce the set 

Iso(g) := {Q G L2((0, l),M2) I sd(g) = sd(g)} 

of potentials isospectral with a given Q and denote by V the set of all 
potentials of the form (13.41) . i.e., 

V := {P = iPij)lj=i I Pij e L2,m(0, 1), pii = 0, pi2 = P2i}- 

The following result is crucial in constructing a pencil Tp ^ with the 
given spectral data. 

Theorem 5.2. Assume that Q E Qh is such that X = is an eigenvalue 
of the Dirac operator ^{Q). Then there exists a unique P E V such 
thatIso{Q)nV = {P}. 

Proof. We divide the proof into three steps. First we prove that there 
is a unique P E V with the property that the transformation opera- 
tor ,^{P,Q) between i{P) and i{Q) on Vq is just the operator ^ of 
multiplication by the matrix- valued function R of (I4.2p - (I4.3I) . On the 
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second step we show that this P is indeed isospectral with Q. Finally, 
we explain why there is no other potentials in lso{Q) fl V. 
Step 1. We write Q G Qh as 

\Q2 -qi 

with real- valued gi and g2 in -^^2(0, 1). Let 6 be an absolutely continuous 
function and R = e^"^; then R commutes with J and one finds that 

R^^ (j— + q)r=J— + R-^JR! + R-^QR. 

\ ax J ax 

Therefore 3?i{P) = £{Q)M for a unique matrix potential P equal to 
(5.3) 

P = R-^JR! + R-^QR 

— (h ni\ T , f Qi cos 29 — q2 sin 29 gi sin 2^ + g2 cos 26* \ 
-[h-t>)I + \^ gisin2^ + g2COs2e -gi cos 2^ + ga sin 2^ J ' 

The potential P defined by (15.31) belongs to the set V if and only if 

(5.4) -^' + gicos2^-g2sin2^ + /i = 0. 

There exists a unique solution of this equation over [0, 1] satisfying 
the initial condition 9{0) = 0. We call this solution 6*2 and determine 
the corresponding potential P & V through (15. 3p with 9 = 92- A 
straightforward calculation shows that this ^2 smd 9i = verify the 
relations (14.31) : in particular, 9^ = h — \p22- By Theorem 14. ![ the 
operator ^ of multiplication by the matrix-valued function R = e^^"^ 
of (14.21) is indeed the transformation operator between i{P) and i{Q) 
on Vq. 

Step 2. Next we claim that ^2(1) = tttt, for some n G Z. By construc- 
tion, M-^i{Q)^ = i{P), so that the operator ^{P) := M'^^{Q)^ is 
a Dirac operator defined via ^{P)u = i{P)u on the domain consisting 
of those u = (mi,M2)* G L2((0, 1),C2) for which ^u G dom^(Q). The 
assumption that A = is an eigenvalue of the Dirac operator ^{Q) and 
the similarity of ^{Q) and &{P) imply that A = is in the spectrum 
of ^(P), and we denote by u° = (mi,M2)* the corresponding eigen- 
function. As -R(O) = / and v° = (fi,f2)* := ^vP is in the null-space 
of ^(Q), we see that ^2(0) = 1-2(0) = 0, and thus M2 = by (Q. 
Therefore, 

.0/1A _ U(^^,.G(^^ _ /^Ml(l)cOS^2(l) 



v"(i) = i?(i)u7r , , , ^ , , 
^ ' ^ ' ^ ' ' -ni(l)sin^2(l) 
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As Ui(l) 7^ and ^2(1) = 0, we conclude that sin6'2(l) = and thus 
that ^2(1) = 7™ for some n G Z. 

Henceforth we have constructed a potential P (z V and a unitary 
operator ^ of the form fl4.2l) - fl4.3p with ^2(1) = 7m for some n G Z 
so that M is the transformation between £(P) and 1{Q) on Vq. By 
Theorem I4.5[ P then belongs to Iso(Q). 

Step 3. Finally, suppose there is another Pi in the set Iso((5) H V. 
By Theorem 14.51 the transformation operator ^(Pi, Q) between ^(Pi) 
and 1{Q) is equal to the operator ^i of multiplication by a matrix- 
valued function e'^^^^^-'^"'"''^*^^)'^, where "di and i?2 are given as in fl4.3p . 
but with P replaced by Pi. Since "(^i = for Hermitian Pi and Q with 
real entries, we conclude that Pi = P by step 1, thus completing the 
proof of Theorem 15. 2[ D 

Theorem 15.21 implies that there is a unique P E V such that 
sd(P) = sd*, for every augmentation sd* of the set sd. As explained 
at the beginning of this section, the Dirac operator ^(P) with every 
such P yields a pencil Tp^q with spectral data sd, thus establishing 
Theorem 12. 7[ 

However, as the augmented set sd* depends on the arbitrary choice 
of tto > 0, different Uq lead to different P eV, and, plausibly, to differ- 
ent quadratic operator pencils Tpg. This uniqueness issue is discussed 
in detail in the next section. 

6. Reconstruction of the pencil: Uniqueness 

In this section, we prove Theorem 12.51 and thus complete the study 
of the inverse spectral problem for quadratic pencils Tp g. Namely, we 
show that the matrix potentials P E V constructed in the previous 
section lead to the same p and q for all choices of the positive parame- 
ter ttQ. 

We again start with an arbitrary element sd G SD, set Aq := 0, 
choose an arbitrary positive number ao, and augment sd with the 
pair (Ao,ao)- Then we construct a shifted AKNS potential Q E Qh 
and the corresponding potential P G Iso(Q) fl V whose spectral data 
coincide with the augmented set, see Proposition 15. II and Theorem 15.21 

Choose now a positive number ag different from a^. Augmenting 
the set sd with (Ao,ao)) we obtain spectral data for another Dirac 
operator ^{Q) with potential Q G Qh, and construct the corresponding 
potential P G Iso(Q) nV. 

It turns out that the potentials Q and Q are related via the so-called 
double commutation transformations, see details in [50], Sec. 3] and [4j. 
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Proposition 6.1. Suppose that two potentials Q and Q from Qh are 
as described above, i.e., that the spectra of the corresponding Dirac 
operators 3i{ff) and S>{Q) coincide and the norming constants an and 
OLn only differ for n = 0. Then 



(6.1) 




Q = Q + Q., 




where 








(6.2) 


Q.- 


= c{x, a^)[v{x)v^ (x) J — Jv(x)v*(x)], 




(6.3) 


c(x, a*) ; 


a* 1 


1 


1 + «*/o^vt(s)v(s)rfs' "* ' ao 


ao 



and V is an eigenfunction of the operator S!{Q) corresponding to the 
eigenvalue Aq = 0. 

Recalling the way the operators &{Q) and ^{P) are related and 
Remark l3.3l on the form of the eigenvector u of &{P) corresponding to 
the eigenvalue Aq = 0, we can write a more explicit formula for Q^,. 

Indeed, by Theorem 14. 5l the transformation operator ^{P, Q) for the 
Dirac differential expressions i'(-P) and 1{Q) on the set Vq is just the 
operator M of multiplication by the matrix-valued function R = e^^'^, 
with 62 being the solution of (15. 4p satisfying ^(0) = 0. Therefore 
V = ^u = e^^"^u; as u = {ui, 0)* and {e^^'^Y = e~^^"^, we easily compute 
that v*v = u*u = uf and 

uu* J — Juu* = u^Ji, 

with 

' 1 



Observing that the matrices J and Ji anticommute, we conclude that 
e^^'^Ji = Jie"^^'^; henceforth, 

vvV - Jvv* = e^2'^[uu* J - Juu*]e-^2'^ 

Set w{x) := 1 + a^ ^ u\{s) ds] then c{x,a^)u\{x) = —w'{x)/w{x) = 
— [\ogw{x)Y, thus resulting in the following form of the potential Q*. 

Corollary 6.2. For the potentials Q and Q of Proposition l6J\ rela- 
tion (16. ip holds with 

Q,{x) = -[\ogw{x)]'e^''^^^'Ji 

(6-4) ., ,,,,/ sin 2^2(3;) cos2^2(x) 

= -log^ a; ' no ) \ ■ oQ( \ 

^ V cos2t/2(x) — sm2fc'2(x) 
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where 62 is the solution of ( 15 ■4p satisfying 6{0) = 0, ui is the first 
component of the eigenvector u of ^{P) corresponding to the eigen- 
value Aq = 0, and 

(6.5) w{x) = 1 + a^: I u\{s) ds. 

Jo 

Now we use the explicit formulae (15.31) and (15. 4p determining the 
potential P from Q and the analogous formula for P and Q to derive 
the crucial result relating P and P. 

Lemma 6.3. For the entries pij and pij of the matrices P and P 
constructed above, the following relations hold: 

P22 = P22, P12 = P12 - (log w)', 
with the function w of ( 16. 5 P . 

Proof. First we recall that 62 is the unique solution of equation ( 15. 4p , 

-6' + qi cos 26 - q2 sin 26 + h = 0, 

satisfying the initial condition 6{0) = 0; here gi and g2 are the entries 
of the AKNS part Qq of the potential Q as in (15. 2p . Likewise, 62 is the 
unique solution of 

-6' + gi cos 2^ -q2sm26 + h = 

satisfying 6(0) = 0, with gi and g2 having similar meaning. Equal- 
ity (16. 4p together with Proposition 16.11 allow to recast the latter equa- 
tion for 62 as 

-6' + gi cos 26 - gs sin 26 + h + (log w)' sm{26 - 262) = 0. 

Observe that 6 = 62 is a solution of this equation satisfying the initial 
condition 6{0) = 0; therefore, uniqueness of solutions yields 62 = 62- 

The potential P is related to Q = Q + Q* through a formula analo- 
gous to ( 15. 3p . i.e., 

P = R-^JR' + R'\Q + Q,)R, 

with R = e^^"^ = e^'^'^ = R. Therefore we find that 

P-P = R-'Q,R = e'''^[~{\ogwye''-^ Jie-''-^]e''^ 

= -(l0gU7)'Ji. 

As a result, P22 = P22 and pu = P12 — (logtf)', and the lemma is 
proved. D 

Corollary 6.4. The potentials P and P constructed above generate the 
same operator pencil Tpg. 
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Proof. In view of (15.1 p and the above lemma it remains to show that 



-Pl2+P'l2 = -P'l2+Pl2- 



By Remark [3.31 we have Ui{x) = exp{J^ pi2{s)ds}, so that u[ = puUi. 
Since w' = a^u\^ w" = 2a^u[ui = 2pi2w', and 



(log^)" = ^ - (^) = 2p,2i\ogwy - [i\ogw)f, 
upon substituting pi2 = pu — (logw)' we find that 

-Pi2 +P12 = -P'i2 + (logw)" + P12 - 2pi2(logw;)' + [(logw)']^ 

= -P'i2 + Pu 
as claimed. The proof is complete. D 

Proof of Theorem \2.5[ Suppose there are two pencils, T = Tp^g and 
T = Tp^^, satisfying assumption (A) and having the same spectral data 
in SD. As explained in Section [3l these pencils lead to two Dirac 
operators S){P) and S^{P) with some potentials P and P in V. 

The spectral data for S){P) and ^{P) can only differ at the norming 
constant for the eigenvalue A = 0; denote these norming constants by 
cto and do respectively. 

Now take do = do and construct the potential P G P as explained at 
the beginning of this section. By Corollary 16. 4[ P and P generate the 
same pencil T. On the other hand, the potentials P and P are isospec- 
tral and belong to V and thus coincide by Theorem 15. 2[ Therefore, T 
and T coincide as well, and the proof is complete. D 

7. Reconstruction algorithm and some extensions 

The proof of existence theorem (Theorem 12. 7p contains explicit steps 
forming reconstruction algorithm. Namely, given an arbitrary ele- 
ment sd of SD, we construct a quadratic pencil Tp^q — i.e., the Sturm- 
Liouville eigenvalue problem (11. II) with potentials p G L2,r(0, 1) and 
q G W2~]^(0, 1) — in the following way: 

(1) fix the enumeration (A„,a„), n G Z*, of the pairs (A, a) in sd 
so that A„ increase, A_i < 0, Ai > 0, and determine the shift h 
from the asymptotic representation of A„; 

(2) augment the set sd with a pair (Ao,q;o), where Aq = and ao 
is an arbitrary positive number; 

(3) construct a Dirac operator Si{Q) with potential Q in the shifted 
AKNS class Qh whose spectral data coincide with the aug- 
mented set sd* (see Proposition 15.10 : 

(4) find the corresponding potential P G Vr\\so{Q) via (15. 30 - 05. 4p : 
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(5) compute the potentials p and q using formulas f IS.ip . 

We finish the paper with several comments. Firstly, the above algo- 
rithm can be used to reconstruct the potentials p and q under different 
assumptions on their regularity. Namely, if p and a primitive r of g 
belong to Ls{0, 1) with s > 1, then the corresponding set SD of spec- 
tral data allows an explicit description (the only difference with the 
s = 2 case is in the decay of the remainders A„ and a„) and the steps 
of reconstruction are as above; cf. the characterization of the spectral 
data for the corresponding class of Dirac operators in |3]. Similar char- 
acterization of the set SD is available if p and r belong to W^2*(0, 1) 
with s > 0; cf. the results of [T71I1B] on eigenvalue asymptotics for 
Sturm-Liouville operators with potentials in Sobolev spaces. 

Secondly, the approach described is not restricted to the Dirichlet 
boundary conditions and can be used to reconstruct energy-dependent 
Sturm-Liouville equations under quite general separated boundary 
conditions. 

Finally, reconstruction from different sets of spectral data (e.g., from 
two spectra, or from Hochstadt-Lieberman mixed data) using this 
method is also possible and will be considered elsewhere, cf. [H]. One 
can also get a Hochstadt-type results [Hj on explicit form of the po- 
tentials when only finitely many spectral data are changed. 
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